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Abstract 

In double-field inflation, which exploits two scalar fields, one of the fields rolls 
slowly during inflation whereas the other field is trapped in a meta-stable vacuum. 

The nucleation rate from the false vacuum to the true one becomes substantial enough 
that triggers a first order phase transition and ends inflation. We revisit the question 
of first order phase transition in an “extended” model of hybrid inflation, realizing the 
double-field inflationary scenario, and correctly identify the parameter space that leads 
to a first order phase transition at the end of inflation. We compute the gravitational 
wave profile which is generated during this first order phase transition. Assuming 
instant reheating, the peak frequency falls in the 1 GHz to 10 GHz frequency band and 
the amplitude varies in the range 10 -11 < Hqw^ 2 IS 10~ 8 , depending on the value of 
the cosmological constant in the false vacuum. For a narrow band of vacuum energies, 
the first order phase transition can happen after the end of inflation via the violation 
of slow-roll, with a peak frequency that varies from 1 THz to 100 THz. For smaller 
values of cosmological constant, even though inflation can end via slow-roll violation, 
the universe gets trapped in a false vacuum whose energy drives a second phase of 
eternal inflation. This range of vacuum energies do not lead to viable inflationary 
models, unless the value of the cosmological constant is compatible with the observed 
value, M ~ 10 -3 eV. 


1 e-mail: amjad.ashoorioon@pliysics.uu.se 



1 Introduction 


The mechanism that ends inflation is still an open question in cosmology. In fact it was 
the achilles heel of the old inflation model pQ, as the universe would never recover from 
inflation. Accelerated expansion of the inflationary false vacuum background pushes the 
walls of true vacuum bubbles exponentially apart such that they never manage to coalesce 
and end inflation. New inflation [2], with the evolution of the inflaton, allowed for another 
mechanism of termination of inflation. Fast-roll evolution of the scalar field at the end of 
inflation can violate the slow-roll condition and end inflation. However there still remained 
to address how the universe has been reheated from the supercooled phase. For that one 
has two assume that the inflaton is coupled to other degrees of freedom and decay of the 
inflaton to these degrees of freedom happens either perturbatively or nonperturbatively [3j. 
For this to happen one often has to fine-tune the couplings or the bare-masses of these new 
degrees of freedom. 

End of inflation through bubble nucleation has this advantage that reheating the universe 
happens from the collision of true vacuum bubble walls and their natural conversion to 
radiation. The idea of new and old inflations were combined in “double-field’’ inflation m 
in which one of the fields rolls during inflation, as in slow-roll inflation, and the second field 
is initially trapped in the meta-stable vacuum (as can be seen in our example, this vacuum 
could be the only existing minimum in the beginning of inflation. The true vacuum could 
develop as inflation proceeds). As the first field rolls, the nucleation rate from the meta¬ 
stable vacuum to the true one becomes large enough that bubbles of true vacuum can indeed 
percolate and end inflation. This is very similar to the false-vacuum dominated Hybrid 
inflation PC! where the tachyonic instability of the waterfall field is replaced with a first 
order phase transition. A prototype of such a potential takes the form 

V(<f>^) = V t 3 + H 1 (0)+H 2 (0,0), (1-1) 

where Vo is the vacuum energy which is constant and 0 is the rolling field whose potential 
V\ (0) along with the vacuum energy drives inflation. 0 is the field which facilitates the first 
order phase transition to the true vacuum. 

In this paper we focus on an example of such double-field inflationary models, realized 
in the context of extended hybrid inflation [7j. Previously, [8j had tried to identify the 
parameter space which allows for exit from inflation with a first order phase transition. In 
their calculation, they had assumed that the pre-exponential factor in the nucleation rate 
per unit four-volume, F = Aexp(— Se), is given by the amplitude of the potential in the 
false vacuum, Vo. However, a better approximation for such a pre-exponential factor is 
given by the mass of the field in the meta-stable vacuum mm- We noticed that, upon 
such a replacement for the pre-exponential factor, the parameter space of models that exit 
inflation through a first order phase transition expands substantially. Having found the new 
parameter space, we calculate the amplitude of the gravitational radiation from such a first 
order phase transition. Like any gravitational wave profile which is generated from a first 
order phase transition, the spectrum has an inverted V-shape which is determined by the 
peak frequency and the amplitude at the peak frequency. The peak frequency also depends 
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on the reheating temperature after the first order phase transition, but if all the energy in the 
false vacuum phase is converted to radiation, the peak frequency falls somewhere in the GHz 
range, which is in the frequency range probed by Chongqing University [13] but below its 
current sensitivity limit by 1CU 2 — 1CT 5 . Nonetheless, it is hoped that future improvements 
of the detectors improve their sensitivity and close the gap between the predicted signal and 
the current sensitivity limit. The shape of the produced gravitational spectrum can be used 
to distinguish the mechanism of termination of inflation from the parametric resonance [3] 
which leads to gravitational wave spectrum with different profile mm- 

We also investigate the gravitational profile from inflationary models that exit from in¬ 
flationary phase through violation of slow-roll but after the termination of inflation the 
nucleation rate from meta-stable vacuum to the true vacuum becomes large enough that 
first order phase transition occurs. If the reheating is instantaneous and efficient, the peak 
frequency from such phase transitions lies in the tera to peta-Herz band and its amplitude is 
much smaller than the models that exit from inflation through first order phase transition. 
At the moment there is no planned probe that aims to the detection of such high frequency 
gravitational waves. 

In this work, we chose the rolling potential to be the quadratic potential for simplicity. 
Also, the large-scale predictions of the model were very close to m 2 0 2 prediction in the ns —r 
plane in the limit of vanishing vacuum energy. This part of the predictions of the model 
was still within the 2 a limit of Planck 2013 results [16]. After the revelation of Planck 2015 
results EH, the prediction of the model now falls outside the 2a region in the n s — r plane 
excluding the running of scalar spectral index from the parameters. Nonetheless, the region 
close to the prediction point of m 2 </> 2 is still within the 3cr region. Allowing for the running 
of scalar spectral index, this region comes back to the 2a C.L. region. It is expected that 
taking the rolling potential to be a low energy scale model, like the hilltop model |T8] or 
the Starobinsky model jT2], the predicted values for r for the models exiting inflation with 
a first order phase transition have a larger overlap with the 2cr region of Planck 2015 data. 

The outline of the paper is as follows. First we will explain the setup that can realize the 
idea of double-field inflation. We identify the parameter space for which inflation can end 
through a first order phase transition and then calculate the power spectrum of gravitational 
waves produced during such phase transitions. Then we identify the region of parameter 
space where phase transition happens after the termination of inflation through slow-roll 
violation. We also compute the gravitational signature of such phase transitions. At the 
end, we show that the universe gets trapped in the meta-stable vacuum if the vacuum 
energy is less than a threshold. We comment on the viability of such inflationary models 
and then conclude our paper. 
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2 Extended Hybrid Inflation 


We focus on the extension of hybrid Inflationary potential mm given by 0 
V (0, 0) = ^ A M 4 + im 1 2 0 2 + 

|A0 4 -^y M^ 3 + ^AW + ^aMV (2.1) 

which is supplemented by the cubic term for the field 0 to provide the possibility of first 
order phase transitions at the end of inflation. For large values of 0, the potential has one 
minimum in both 0 and 0 direction. As inflation proceeds and 0 rolls toward its vacuum, 
second minimum along the 0 direction develops if y 2 > 4a , y at 

,9 „ , 9 y 2 — 4a A 

inflection = M 2 ^ (2.2) 

The two minima are separated by a barrier that allows for a first order phase transition from 
the meta-stable vacuum to the true one if the nucleation rate is substantial and when it is 
energetically favorable. Mapping the potential to (ED, Vi( 0 ) is the vanilla chaotic model 
m 2 0 2 and the expression in the second line plays the role of 10 ( 0 , 0 ), which couples the 
inflaton to the field that creates the false vacuum in the 0 direction and facilitates the phase 
transition. In principle, 10(0) determines the predictions of the model at large scales in the 
limiting case that the cosmological constant goes to zero and can be chosen such that the 
model is compatible with the CMB observables at cosmological scales [16]. 

The probability of phase transition is given by the nucleation rate divided by the 4- 
dimensional Hubble volume 

P=|I- ( 2 - 3 ) 

where T is the nucleation rate [ 120 ] 


T = A exp(—S e) (2.4) 

where S'e is the Euclidean four-dimensional action for the solution that interpolates between 
the two minima. For a first order phase transition with a quartic polynomial potential, the 
numerical results were fit by [ 21 ] to have the form 

47T 2 

Se = 3 ^-(2 - 5)- 3 (ai5 + a 2 d 2 + a 3 S 3 ), (2.5) 

where a\ = 13.832, a 2 = —10.819, a 3 = 2.0765, and S is a function of 0 2 , 

_ 9Act 9AA'0 2 
' “ + • 

1 For the status of hybrid inflation with quadratic rolling potential after Planck 2013 data release, please 

see [2]. In particular the paper discusses an interesting scenario in which the first part of the required number 

of e-folds, needed to solve the problems of Big Bang, is provided from the chaotic phase and the second part 
is resulted from the vacuum dominated phase of hybrid inflation. The model can render lower scalar spectral 
indices but larger amount of tensor-to-scalar ratio in comparison with the purely quadratic potential. 
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The allowed range has 0 < 5 < 2. Prefactor A has dimension mass 4 . In |8], this prefactor is 


taken to be equal to the cosmological constant in the false vacnnm, 


AM 4 


However as pointed 


out by Linde in [ID], if M 2 > 2H 2 , the pre-exponential factor A is closer to M 4 . As we will 
see, upon this change, the parameter space in which exit from inflation can happen through 
a first order phase transition expands substantial!}!!. 


For illustration we focus on the following set of parameters 


a = 0.01, A = 1, & A' = l. (2.7) 

All massive parameters are in unit of Planck mass, mp = G -1 / 2 = 1.2209 x 10 19 GeV. For 
the potential (12.ip to have a zero vacuum constant after the phase transition, the following 
relation between the parameters should hold 


2 -a 3 + 9ckA 2 + ( a 2 + 3A 2 ) 3/2 

7 =- y -— ( 2 - 8 ) 

As in [8], we take the required number of e-folds to solve the standard big bang cosmology 
problems to be 55. Inflation can also end through end of slow-roll which is parameterized by 


_ m.pi /H/A 2 m 4 0 2 mp 

167t\Vi/ 7t(AM 4 + 2?n 2 0 2 ) 2 
which has the biggest solution 

2 m 2 mp 1 + mmp y / m 2 m| 1 — 8tiXM 4 — 47tAM 4 

e 87 t m 2 


(2.9) 


( 2 . 10 ) 


as the physically viable solution. 

In order for inflation to end through phase transition rather than slow-roll violation, 
0 pt > 4>e- The number of e-foldings could be obtained as a function of the inflaton field. As 
stated above, we take the required number of e-folds equal to 55: 


-W*,0 P t) 


87T 

rrip 



V 1 


Vi 


7 d(f> 


2vrA- 


M 4 


rnrni 


■ In ■ 


+ 


2vr 


mi 





( 2 . 11 ) 


From the above one can obtain 0*. The scalar spectral index and the tensor over scalar ratio 
then could be obtained at 0* using 


ns — 1 = —6e + 2r ]; 

r = 16e, 


where e is given above and 


Trip V" m 2 mp 

87t V 2tt(XM 4 + 2m 2 0 2 ) 


( 2 . 12 ) 

(2.13) 


(2.14) 
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Figure 1: The behavior of m (left) and n$ (right) as a function of M in the region of 
parameter space that inflation ends with a first order phase transition. 

3 End of Inflation with a first Order Phase Transition 

In this section we focus on the region of parameter space where before the violation of the 
slow-roll parameter, the nucleation rate becomes substantial enough that leads to percolation 
of true vacuum bubbles in the sea of inflating false vacuum. We assume that the pre¬ 
exponential factor, A is of order , as in [TO], which makes our calculations distinct from 
the analysis of [Sj. As we will see, this assumption will expand the region of parameter 
space in which end of inflation happens through a first order phase transition. This part 
of parameter space is part of the region in which the vacuum energy contribution to the 
potential is comparable to the energy density of the inflaton field </>. Detailed computation 
shows that only in the region where 

4.97 x KT 4 < M < 2.66 x 10" 3 , (3.1) 

first order phase transition precedes the slow-roll violation. For M > 2.66 x 10 -6 , there is 
no solution to the density perturbation amplitude normalization and hence no viable model. 
One can match the amplitude of density perturbations with the COBE normalization which 
fixes the mass parameter of the scalar field, m. In the left plot of figH] we have graphed how 
the mass parameter varies as a function of M . The obtained range of m is 

1.24 x 1(T 6 < m < 6.33 x KT 6 . (3.2) 

Within the range (13.ip . ns and r vary in the following ranges 

0-91 <n s < 0.97, (3.3) 

0.15 <r <0.55, (3.4) 

as in the right graph in figH] and left plot of figM We have also plotted the predictions of 
this region of parameter space in the n$ — r plane. As can be seen from the plots, with 
increasing M, initially n$ decreases and r increases. However there are turning points: 

2 If M 2 < 2H 2 , the prefactor is of order H 4 [IT]. 
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Figure 2: The behavior of r vs. M (left) and vs. ns (right) in the region of parameter space 
that inflation ends with a first order phase transition. 


around M ~ 2.2 x 10 -3 , ns starts to increase and around M ~ 2.4 x 10~ 3 , r starts to 
decrease. The loci of predictions of the model is designated in the ns — r plane in the right 
plot of figI2J For small values of M, the predictions of the models for large scale fluctuations 
is very close to the chaotic quadratic potential. 

With the choice of parameters as in (12.71) . the true minimum in the 0 direction appears 
very late, i.e. towards the end of inflation. There is no minimum for 0 that a bubble of true 
vacuum forms and lead to open inflationary scenario, as in [12] , 


If the energy stored in the inflaton potential is completely transformed to the radiation, 
i.e. reheating is instantaneous, one can calculate the reheating temperature at the end of 
inflation 


T* = 


90 H) 
87 r 3 g* 



(3.5) 


where g * is the the total number of relativistic degrees of freedom at temperature T* 
which we take to be g * ~ 106. We have plotted T* vs. M in the left plot of fig]3] If reheating 
is not efficient, T* is smaller than what is given in (13.51) . This temperature determines the 
peak frequency of the gravitational wave spectrum generated through the first order phase 
transition. If one assumes instant reheating, the reheating temperature varies almost linearly 
in the range 10 15 GeV and 10 16 GeV vs. M in the range (13.ip . 


3.1 Gravitational Wave Spectrum from Models with First Order Phase Tran¬ 
sition at the end of Inflation 

Since inflation supercools the universe, one can exploit the formalism of gravitational wave 
generation from first order phase transition at zero temperature. The numerical computa¬ 
tions for bubble collisions from a first order phase transition was initially done by [24j for two 
bubbles and generalized for more bubbles in [25]. The spectrum has the shape of asymmetric 
of inverted V around a peak frequency, / m , which decays like f 28 and / _1 respectively at 
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Figure 3: The behavior of T* and (■Jy) vs. M in the region of parameter space that 
inflation end with a first order phase transition. 


smaller and larger frequencies. The peak frequency f m today, after the post-inflationary 
redshifting is taken into account, is given by [2TJ 


/. 


3 x 10” 10 



(iiv) 


and the amplitude at this peak frequency is given by 



^gw^ 2 (/?7 


10“ 6 




(3.6) 


(3.7) 


where 0 is 


gLSe dS-E dcf) 
dt d<f> dt 


(3.8) 


0 _1 is a measure of how fast the first order phase transition takes to complete. For our 
computations to be reliable, we expect this time to be much smaller than the expansion rate 
of the universe, i.e. fd/Hf 1, where 77/ is the Hubble parameter at the time of phase 
transition which coincides with the end of inflation. ^ is the velocity of the scalar field, 0, 
which can be found during inflation through the following relation [26] 


2 2e(0)W(0) 

3 - e(0) 


(3.9) 


where e is the first slow-roll parameter. We have plotted log vs. T7 in the range (13. ip 

of M that leads to exit from inflation with a first order phase transition. As can be seen 
from the right plot of fig J3J for smaller values of M in the range, phase transition is quite fast 
in comparison with the expansion rate of the universe, 0/77/ ~ few x 100. As M increases, 
phase transition takes longer to complete in comparison with the expansion time of the 
universe. For the maximum value of M in the range, M = 2.66 x 10~ 3 , 0/77/ = 6.11, which 
is nonetheless fast enough to validate our computations. The intensity of the gravitational 
waves at the peak frequency is only dependent on this parameter, 0/77/. The faster the phase 
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Figure 4: The peak frequency of the gravitational wave spectrum and its amplitude as a 
function of M for models that exit inflation with a first order phase transition. 


transition, the smaller the amplitude of the produced gravitational waves. We have plotted 
log(f2cw^ 2 ) vs. M which increases in the range (13.11) from 1.07 x 10~ n to 2.63 x 1CT 8 , see 
the right plot in figldj However the peak frequency, besides (3/Hf depends on the reheating 
temperature too. Assuming that the reheating is instantaneous and the total energy density 
of potential at the end of inflation transforms to radiation, we have plotted the peak frequency 
as a function of M in the range (13.11) . see left plot of fig [I] For the smallest value of M in the 
range, M ~ 5 x ICG 4 , peak frequency, / m is around 1.63 x 10 10 Hz whereas for the largest 
value of M, / m ~ 1.75 x 10 9 Hz. Such a gravitational wave profiles lies in the frequency band 
of Chongqing High Freqeuncy Gravitational (HFGW) probe but below its current sensitivity 
limit by a factor of 10~ 2 — ICC 5 [13] . It is hoped that the improvement of the detector in 
future close the gap between the expected signals and the current sensitivity limit. 

If the reheating coming from bubble wall collision is not instantaneous and a phase of 
non-radiation domination intervenes the end of first order phase transition and the radiation 
domination, T* will be lowered. If this efficiency factor is assumed to be y, where y < 1, 
the amplitude at the peak frequency will be lowered by a factor of y 4 [30]. For example if 
y ~ 0.01 — 0.1, the frequency range will be shifted within the sensitivity band of Birmingham 
HFGW probe [27]. If the reheating temperature is a factor of y ~ ICG' — ICG 10 smaller, the 
peak frequency will shift to the sensitivity band of DECIGO [28] and BBO [29] . However in 
all these cases, the amplitude of the signal will be lowered such that the signal could not be 
observed by any of these probes. 

With the choice of V\ ((j)) as the quadratic potential, most of the parameter space of 
inflationary models that exit from inflation with a first order phase transition is ruled out. 
There is a small region of parameter space which has predictions very close to m 2 </> 2 . This 
region in particular is of interest for the CMB polarization probes like the BICEP2 [22] or 
the future ones like the CMBPol [23]. Since these models exit inflation through first order 
phase transition which is accompanied by bubble collision, they can leave an extra signature 
of gravity waves at higher frequency scales. We expect that the predictions of the model for 
r could be lowered if one would replace V\ with another model, like the hilltop model, [IB], 
or models with smaller energy scales that naturally predict a lower value for r. 
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Figure 5: m (left) and ns (right) vs. M for the models for which first order phase transition 
happens after the end of inflation. 

4 First Order Phase Transition after the End of Slow-roll Inflation 


It is possible that even though inflation ends through violation of the slow-roll inflation, the 
first order phase transition happens after the end of inflation, when <f> becomes closer to the 
meta-stable vacuum, 0 = 0. In fact if 

3 x 10“ 6 < M < 4.88 x KT 3 4 (4.1) 

phase transition happens after the end of slow-roll inflation. The lower bound of the above 
interval is intriguingly very close to the scalar field mass in the rolling direction. To be able to 
determine this range of M, we solved the equations of motion and computed the evolution 
of the scalar field after inflation numerical!}!!. As before, we have to pick up only the 0 
solution to the nucleation rate equation that 0 < 5 < 2. We have plotted the variation of m 
as a function of M in the range (14.2ft . see the left plot in fig® Assuming that the reheating 
is instantaneous, we have also plotted the reheating temperature in this range of M, see the 
right plot in fig® In both cases, parameters m and T* are almost constant for smaller values 
of M in the range, but they gradually rise as M increases. As expected, contrary to the 
behavior of r, ns decreases with the rise of M, please see fig. ®and® We have also plotted 
the behavior of ns vs. r as M increases, see fig® The behavior of log vs. Af, where 
Hf is the Hubble parameter when phase transition completes can be seen in right plot in fig. 
® As can be seen from the plot, with the increase of M, phase transition becomes slower. 
Nonetheless, in general, phase transitions happen much faster after inflation in comparison 
with the inflationary models that end via first order phase transition. This is partly due to 
the fact that the Hubble parameter after inflation ends, is smaller than its corresponding 
value at the end of inflation. As M enhances in the interval, (14.20 . (3/Hf decreases from 
4.21 x 10 6 to 90321.6. We have checked that in this range of M, the effective mass of 0 field, 

3 In the range 

4.88 x 1CT 4 <M< 4.97 x 1CT 4 , (4.2) 

the constraints on the number of e-folds, N e = 55, with end of inflation given by the slow-roll violation did 
not yield a real solution for to. One can satisfy the constraint equations with less number of e-folds though. 
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Figure 6: r vs. M (left) and ns (right) for the models for which first order phase transition 
happens after the end of inflation. 


M 2 > 2 H 2 and thus phase transition happens via Coleman-de Luccia instanton transitions 

m- 

Using Eqs. (j3.6[) and 03.7H . we have calculated the peak frequency of the gravitational 
wave spectrum and its intensity, Hcw^ 2 , as a function of M, please see the plots in figJHl As 
M increases, f m varies from 1.8 x 10 14 to 4.95 x 10 12 Hz, which is well outside the frequency 
band of any currently planned probe. Of course this frequency range is obtained assuming 
that the reheating is instantaneous and efficient. If the reheating temperature is smaller than 
its instantaneous value by a factor of x — 10~ 3 — 10~ 4 , the peak frequency range is moved 
to the sensitivity band of Chongqing high frequency gravitational wave probe. However, 
the intensity will then gets suppressed further by a factor of 10 -12 — 1CU 16 which makes the 
signal too small to be detected. 




Figure 7: T* (left) and log (/3Hf) (right) vs. M for the models for which Erst order phase 
transition happens after the end of inflation. 
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Figure 8: The peak frequency (left) and the amplitude at the peak frequency (right) vs. M 
for the models of inflation in which first order phase transition happens after inflation. 


5 Trapping in the Meta-Stable Vacuum and Eternal Inflation 


For smaller values of M, M < 3 x 10 -6 , the nucleation rate as given by Coleman-De Luccia 
transition rate, is too small to allow for the first order phase transition to complete as the 
scalar field, 0, passes through the minimum of the potential, V\. Gradually as the effective 
mass of the -0 field decreases while the 0 field approaches the metastable vacuum, there is 
a chance that M, 2 < 2 H 2 . In particular this can happen for smaller values of a and M. In 
this situation, Hawking-Moss phase transition [3T] appears where an inhomogeneous bubbles 
whose radius is greater than de-Sitter space radius, H^ 1 , will form in the flat space-time. 
The tunneling probability per unite four volume is of order [3TJ 

r = TT 4 exp(-H), (5.1) 


where 


E(0,0) H(O,-0 max )J 

is the difference between the combined gravitational and scalar field action of the -0 = 0 max , 
which is another homogeneous solution apart from -0 = 0. V(0, v k max ) is the local maximum 
of V on a four-sphere of radius, i/“ ax , where H^ ax = 87rH(O,0 max )/3. As 0 = 0 max is 
unstable, the scalar field runs downhill from 0 max , to the global minimum 0 min , afterwards. 
For the potential (12.11) . one can calculate the exponent B and show that 



B — 


3A 2 M 4 


(5.3) 


for a <C 1. Therefore for 


M < 


a 3 / 4 
3 1 / 4 A 1 / 2 


(5.4) 


the Hawking-Moss phase transition rate from the meta-stable vacuum to the the true vacuum 
is very small. In our example where we took, A = 1 and a = 0.01, for M < 6.08 x 10~ 5 such 
phase transition takes a lot of time to complete and basically leads to a self-reproduction 
regime like old inflation [[P]. After the Hawking-Moss transition completes, the 0 field will 
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be homogeneous on scalar of order H 0 , with H 0 
stochastically distributed with the dispersion HH 


\J Nonetheless, it will be 


= ( 0 2 ) 


3g 0 4 

87 r 2 m 4 ’ 


(5.5) 


which in our example is about ~ 10 5 and thus too small to drive another phase of inflation. 

In [ 8 ], the authors claim that they have been able to find a branch that corresponds 
to the vacuum-dominated regime of hybrid Inflation [7] in the Extended Hybrid inflation, 
potential (12.11) . However in hybrid inflation, this regime is obtained assuming that the 
tachyonic instability in the waterfall field ends inflation. In this case, however the waterfield 
mass squared, aM 2 + Xijj 2 is always positive and never becomes tachyonic. For such energy 
scales in the vacuum dominated regime, the slow-roll can never get violated as the first 
slow-roll parameter decreases as inflation progresses. We also showed that Hawking-Moss 
phase transition for such small vacuum energies is not substantial enough to end inflation. 
Therefore, it is not possible to realize the vacuum energy dominated regime of hybrid inflation 
in such extended models. The scalar field gets trapped in the metastable vacuum with no 
graceful exit and therefore this region of parameter space does not yield a viable inflationary 
model. On the other hand, from the phenomenological perspective, if M is of order ~ lO -30 , 
this vacuum energy can be responsible for the current acceleration of the universe. 


6 Conclusion 

Terminating inflation with a first order phase transition has the benefit of reheating the 
universe from the supercooling phase of inflation through the collisions of bubbles of true 
vacuum without invoking and fine-tuning of the couplings of the inflaton to the other fields. 
One can achieve this scenario, modifying the old inflation scenario, with time-dependent 
nucleation rate which is small in the beginning and becomes substantial at the end of infla¬ 
tion. One specific realization of this scenario is extended inflation [32], in which the gravity 
sector of the theory is modified to Jordan-Brans-Dicke theory. Another way of achieving 
this scenario, is having two scalar fields, where one of the fields rolls and the other one is 
trapped in a meta-stable vacuum mu- As the rolling field evolves, the nucleation rate at 
the false vacuum becomes large enough that the condition for percolation of true vacuum 
bubbles holds and inflation ends. We re-examined a model of extended hybrid inflation which 
provides us with such a setup. We noticed that the pre-exponential factor in the nucleation 
rate plays a crucial role in correctly determining the parameter space that allows for a first 
order phase transition in the model. For models that exit inflation with a first order phase 
transition, we computed the peak frequency and the amplitude at the peak frequency, which 
are respectively in the ranges f m ~ 10 9 — 10 10 Hz and Ogw ^ 2 — 10 ” 11 — 10 -8 . The signature 
is in the frequency range of Chongqing HFGW probe but smaller than the current sensitivity 
limit by a factor of 10 ~ 2 — 10 -5 . It is expected that future improvements of the detectors of 
the probe can close this gap between the signal predicted by the model and the sensitivity 
limit [IB] . One should note that the predicted amplitude for the high frequency gravitational 
signal is the same as the gravitational signal generated by parametric resonance in the GUT 
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scale models of inflation [33] [MJ , although the two spectra could be distinguished from each 
other by their shape. In particular the gravitational wave profile generated by the parametric 
resonance has a double-hump feature instead of the single peak profile of the gravitational 
wave signal generated by the bubble collision. 

We also showed that it is possible for some models to make a first order phase transition to 
the true vacuum after the termination of inflation with violation of the slow-roll parameter. 
However the signal of such a phase transition is weak and out of the frequency band of future 
probes. 

The large-scale predictions of the model are very dependent on the potential of the rolling 
field. Here with the choice of quadratic potential for the rolling field, the predictions of the 
model approaches the predictions of m 2 (j) 2 in the ns — r plane, which was still within the 2cr 
limit of Planck 2013 results [16] . As this model produces gravitational waves at the CMB 
scales, the analyzed model produces an observable gravitational waves at both small and 
high frequency range of the spectrum: the small frequency gravitational waves have quantum 
origin but the high frequency gravitational waves have classical origin. After the exposure 
of Planck 2015 results mi, major prediction of the model is now outside the 2cr region in 
the n s — r plane excluding the running of scalar spectral index from the parameters. The 
region close to the prediction point of m 2 0 2 is still within the 3cr region though. If one allows 
for the running of scalar spectral index, this region comes back to the 2 a confidence region. 
One should also note that the predictions of the model at large scales is very sensitive to the 
initial condition for fluctuations. As it was shown in [35], choosing excited initial condition 
for cosmological perturbations generally tend to suppress the tensor-to-scalar ratio. 

ft is expected that if one takes the rolling field potential to be a lower energy scale 
model, like the hilltop [18] or the Starobinsky [19] model, the predicted values for r for the 
models exiting inflation with a first order phase transition have a larger overlap with the 2cr 
confidence region of Planck 2015 results. This is something that 1 will postpone to a future 
publication. 

As the peak frequency of the gravitational spectrum, produced from the first order phase 
transition, depends very much on the Hubble parameter at the time of phase transition, it is 
also interesting to investigate inflationary models that exit inflation with smaller Hubble pa¬ 
rameter and, hence, peak frequency that falls within the sensitivity bands of BBO, DECIGO 
or even Advanced L1GO. 
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